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Abstract.  In  an  earlier  note  Cimbala  and  Park  (1990)  presented  a 
method  for  calibrating  constant  temperature  hot  wires  in  incom- 
pressib'o  gas  flow  which  accounts  for  ambient  temperature  changes. 
Their  method  involves  the  retention  of  the  temperature  dependence 
in  a  form  of  King's  Law.  However,  in  practice  the  method  proves 
diflleult  to  implement  due  to  a  numerical  difficulty  in  finding  a  best 
fit  to  the  assumed  functional  form.  This  note  presents  a  modified 
version  of  the  method  of  Cimbala  and  Park  which  docs  not  seem  to 
have  the  numerical  difficulty,  as  well  as  some  further  comments  on 
the  method. 


1  Introduction 

In  an  earlier  note,  Cimbala  and  Park  (1990)  presented  a 
version  of  King’s  law  which  retains  a  dependence  on  the 
ambient  temperature,  and  they  showed  that  this  method 
allows  reasonably  accurate  computation  of  the  velocity  as  a 
function  of  hot-wire  voltage  over  a  wide  range  of  tempera¬ 
tures.  In  common  with  most  practical  use  of  King’s  law,  the 
calibration  coefficients  computed  using  their  method  are 
dimensional,  and  this  tics  into  a  numerical  difficulty  with  the 
method,  which  is  that  finding  the  values  of  the  calibration 
coefficients  which  give  the  best  fit  to  the  assumed  form  of  the 
calibration  curve  is  numerically  unstable,  and  requires  sub¬ 
stantial  cfTort  in  order  to  obtain  convergence.  By  using  a 
nondimensiunai  form  of  King’s  law,  the  calculation  of  the 
best  fit  is  simplified  by  the  fact  that  the  nondimensional 
coefficients  and  variables  arc  all  of  order  1. 

The  starting  point  is  the  following  form  of  King's  Law: 

1V»  =  fl  +  h  Ke"  (1) 


where  Nu  is  the  Nusselt  number  and  Re  is  the  Reynolds 
number  based  on  wire  diameter  and  the  local  flow  velocity 
U  normal  to  the  wire.  Following  Black  welder  (1981)  and 
modeling  the  hot  wire  as  a  circular  cylinder,  the  Nusselt 
number  is  defined  by 

nl^k  (T„,  -  TJ 

where  q  is  convective  heat  lo.ss,  /„,  is  the  hot-wire  length. 
k  -  k{T)  is  the  thermal  conductivity  and  and  7^  arc  the 
wire  operating  temperature  and  the  ambient  temperature, 
respectively.  Assuming  that  the  heat  loss  is  equal  to  the 
power  dissipated  in  the  wire. 


where  t'„.  is  the  voltage  across  the  hot  wire  and  R„.  is  the 
(nominally  fixed)  hot-wire  resistance.  This  gives 

^  k  (7^  7^) 

As  in  Cimbala  and  Park,  k  and  the  kinematic  viscosity  v  are 
evaluated  at  r=i(T„, -F  TJ  (Sutherland’s  law)  and  at  the 
ambient  pressure.  Substituting  Eq.  2  into  Eq.  1  provides  the 
desired  explicit  dependence  on  the  ambient  temperature. 

Cimbala  and  Park  rewrite  this  equation  in  dimensional 
form  with  5  calibration  parameters,  corresponding  to  a  hot¬ 
wire  offset  voltage,  the  wire  temperature  7;.,,  the  exponent  m. 
and  dimensional  forms  of  the  coefficients  a  and  h.  ’The  offset 
voltage  is  not  a  calibration  parameter  in  the  same  sense  as 
the  other  parameters  as  it  can  be  independently  measured 
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and  is  not  a  property  of  the  hot  wire  itself,  and  therefore  is 
not  used  as  a  parameter  in  the  present  method.' 

If  the  equation  is  evaluated  directly  in  the  form  given 
here,  Nu  and  Re  arc  of  order  1  and  minimization  of 

=  Z  (Nu{Ui,  T„,,  P,;  a,  h,  n,  TJ  -  Nu,f 

is  straightforward,  giving  the  values  of  the  fit  parameters 
(a,  b,  n,  TJ  for  use  in  calibrating  the  hot  wire.  Although  T„. 
is  a  dimensional  value  of  order  100  K,  and  is  used  as  one  of 
the  fit  parameters,  this  does  not  affect  the  minimization 
scheme  since  it  is  combined  into  Nu  which  is  of  order  1.  The 
expression  for  Nu  gives  e„{U,  TJ  and  is  directly  invertible, 

‘  It  is  straightforward  to  include  an  olTsct  voltage  by  replacing  the 
hot-wire  voltage  with  c—e„„,  where  e  and  e^f  arc  the  measured 
hot-wire  signal  and  offset  voltages  respectively,  and  then  minimiz¬ 
ing  with  respect  to  the  additional  parameter  e„„ 


so  U  (e„,  TJ  is  easily  calculated  once  the  values  of  the  fit 
parameters  are  known. 

A  large  quantity  of  data  was  taken  in  air  with  a  commer¬ 
cially  available  tungsten  hot  wire  at  different  ambient  tem¬ 
peratures  spanning  a  range  from  18-30°C,  and  values  of  the 
fit  parameters  were  determined  using  a  simplex  method 
(Press  ct  al.,  1988)  to  search  for  a  multi-dimensional  mini¬ 
mum  of  in  the  (a,  h,  n,  TJ  parameter  space.  The  physical 
properties  of  air  were  computed  by  interpolation  from  the 
tabulation  in  Hilsenrath  (1955).  The  simplex  method  re¬ 
quires  input  of  an  initial  guess  of  several  points  in  the  coef¬ 
ficient  space,  and  a  tolerance  defining  how  close  to  the  actual 
minimum  the  result  must  be.  In  the  present  case,  a  single 
initial  guess  is  provided,  and  the  other  points  are  produced 
by  moving  15%  of  the  initial  value  along  each  axis  in  the 
parameter  space.  The  tolerance  is  defined  as  the  maximum 
allowable  (Euclidean)  distance  from  a  minimum  in  parame- 
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Tabic  I.  Hcsi  fit  parameters  for  the  hot-wire  calibration  data  shown 
in  I  ig.  la 


I'aramclcr  Ucst  fit  value 

«  0.2.1  IX 

!'  0..1I7I 

n  0.4225 

.1X1.0  k"* 


ter  .space,  and  110  '■  was  round  to  give  accurate  results. 
Tlie  fit  parameters  found  using  the  experimental  data  arc 
listed  in  Table  I.  The  value  of  T„.  in  Titblc  I  which  minimizes 
■/'  should  reflect  the  actual  hot-wire  temperature,  which  is 
selected  to  be  500  K  by  a  setting  in  the  hot-wire  circuit,  and 
it  is  uncletir  why  the  value  found  by  the  calculation  is  so 
much  smaller.  The  method  w'as  found  to  be  very  tolerant  of 
the  choice  of  the  initial  guess,  and  for  the  example  data 
always  converged  to  the  same  values  of  the  parameters. 

Figure  la  is  a  plot  of  the  data  in  the  form  of  Ufn  vs. 
^nivii-un.i  •  "dll  evaluated  using  the  coefficients  in  Table  I. 
and  Fig.  lb  is  a  plot  showing  the  relative  deviation  of  Uf^ 
from  •  Tlic  uvcragc  deviation  of  the  data  from  the 

measured  values  was  0.06  m/s.  which  is  similar  to  that  found 
by  Cimbala  and  Park,  and  it  can  be  seen  that  the  relative 
deviation  increases  for  low  velocities.  Further  calibration 
data  for  the  same  hot  wire  was  collected  on  a  different  day 
from  the  data  shown  in  Fig.  la  and  1b,  and  is  shown  in 
Fig.  Ic.  The  result  of  processing  the  data  using  Eq.  1  is  a 
straight  line,  which  may  then  be  corrected.  After  this  second, 
linear  correction  is  applied  the  relative  deviations  (not 
shown)  are  all  le.ss  than  1%  for  U>2  m/s,  with  an  average 
deviation  of  0.04  m/s.  This  demonstrates  that  Eq.  1  is  a  good 
linearizer  of  the  data,  but  it  cannot  be  used  for  velocity 
measurement  without  a  minimal  recalibration  to  determine 
the  linear  correction,  which  is  probably  caused  by  hot-wire 
drift  from  the  effects  of  oxidation  or  other  factors.  A  linear 
correction  was  also  used  by  Cimbala  and  Park. 

The  values  of  the  coefficients  in  Table  1  were  also  used 
to  linearize  data  taken  for  a  different  hot  wire  of  the  .same 
make  and  model  and  Fig.  Id  is  an  example  of  the  fit  for  such 
a  hot  wire.  Again,  after  applying  a  linear  correction,  the 
deviations  in  this  case  are  all  less  than  1%  for  U>2m/s. 
with  an  average  deviation  of  0.05  m/s. 

From  an  examination  of  the  figures,  it  appears  that  there 
is  a  small  but  consistent  curvature  in  the  “linearized"  data 
sets.  In  an  attempt  to  correct  for  this  a  quadratic  term  was 
added  to  the  King's  Law  expression,  resulting  in  the  form: 

i\’ii  =  (i  +  h  Re"  +  (■  Re-" 


The  additional  term  can  be  thought  of  as  the  next  term  in  a 
power  series  for  Nu  {Re"),  and  so  is  not  entirely  without 
foundation.  This  expression  is  also  invertible  in  closed  form. 
Using  the  data  set  of  Fig.  la,  and  the  simplex  method  de- 
■scribed  above,  the  value  of  c  was  found  to  be  5  orders  of 
magnitude  less  than  a  and  b.  and  the  deviation  was  not 
reduced. 

The  experiments  for  which  this  effort  were  undertaken 
have  been  moved  to  a  facility  where  the  temperature  is  typ¬ 
ically  constant  to  within  0.1  K  per  hour,  and  under  these 
circumstances  it  is  simpler  to  compute  (u,  h,  n,  TJ  using  cal¬ 
ibration  data  taken  for  each  experiment,  rather  than  at¬ 
tempting  to  use  a  single  set  of  coefficients  obtained  from  a 
large  quantity  of  data. 

In  conclusion,  the  temperature  dependent  King’s  law 
given  in  Eqs.  1  and  2.  combined  with  a  linear  correction, 
provides  an  effective  means  for  making  hot-wire  measure¬ 
ments  in  facilities  where  variations  in  ambient  temperature 
make  the  more  traditional  method  impractical  or  impossible 
to  use.  The  modified  version  given  in  this  note  is  useful 
because  of  the  ease  with  which  the  fit  parameters  may  be 
computed  and  is  philosophically  satisfying  because  the  coef¬ 
ficients  appearing  in  the  equation  arc  dimensionless.  Treat¬ 
ing  King's  law  as  a  truneated  power  series  for  Nu  {Re")  and 
adding  the  next  term  in  the  scries  was  shown  to  have  no 
significant  effect  on  the  quality  of  the  fit  to  the  data.  Further, 
the  values  computed  for  one  hot  wire  have  been  shown  to 
satisfactorily  linearize  calibration  data  from  a  different  hot 
wire  of  similar  construction,  so  that  extensive  calibration 
measurements  over  a  wide  temperature  range  need  be  per¬ 
formed  for  only  one  hot  wire,  and  the  resulting  calibration 
parameters  then  used  for  several. 
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1  Introduction  \ 

I'xperimcnlal  studies  conducted  in  the  UDIVs  hi’g^fO  to  show 
liitil  the  development  of  plane  mixing  layers  was  iimuenccd 
ni>l  only  by  the  formtition  and  interaction  of  lar^-scalc 
spanvvise  vortices,  but  also  by  a  secondary  structure  ib  the 
form  ofstreamwise  vortices  (Brown  and  Roshko  1974),  ^v- 
er;d  studies  were  subsequently  conducted  on  the  origin  and 
evolution  of  this  secondary  structure,  but  almost  all  of  thenrV, 
were  limited  to  flow  visualization  data  obtained  mainly  in 
the  near-field  region  of  the  mixing  layer  (see  Bell  and  Mehta 
l'W2,  for  a  review). 

The  first  detailed  quantitative  investigation  of  these 
streamwise  structures  was  conducted  in  a  two-stream  mixing 
layer  with  a  velocity  ratio  of  0.6  and  a  Reynolds  number  of 
about  ~2.9x  10'*  (Bell  and  Mehta  1992).  Measurements 
the  mean  streamwise  vorticity  indicated  that  small  distin'- 
banccs  (naturally  present)  in  the  flow  were  initially  ampliflcd 
just  downstream  of  the  first  spanwise  roll-up,  leading  t/ the 
formation  of  streamwise  vortices,  in  agreement  with  /arlicr 
observations.  The  streamwise  vortices,  with  averag^ircula- 
tion  equivalent  to  about  10%  of  the  initial  spanwirc  circula¬ 
tion,  first  appeared  in  clusters  containing  vorti^s  of  both 
signs,  but  re-organized  further  downstream  to  f^rm  counter¬ 
rotating  pairs.  This  vortex  structure  was  foimd  to  grow  in 
size,  .scaling  approximately  with  the  mixinar  layer  vorticity 
thickness,  and  weaken,  the  maximum  mctyl  vorticity  diffus¬ 
ing  as  approximately  l/A"'  *.  The  data ^ggested  that  the 
streamwise  structures  persisted  ihrou^  to  the  far-field  re¬ 
gion,  although  they  were  weak  cnougfyby  this  point  that  the 
mixing  layer  may  be  considered  to  lynominally  two-dimen¬ 
sional.  / 

Mixing  layers  are  known  to^e  very  sensitive  to  initial 
conditions,  and  so  one  questionyWhich  naturally  aro.se  in  the 
above  study  was  to  what  exttmt  the  details  of  the  observed 
three-dimensionality  were  ftufility  dependent,  i.e.  would  the 
nature  of  the  threc-dimensmnality  be  the  same  in  all  wind 
tunnels?  Hence,  the  objective  of  the  present  study  was  to 
establish  the  sensitivity  m  the  three-dimensionality  to  small 

*  Present  address;  School  of  Mechanical  F.nginecring,  Purdue  Uni¬ 
versity,  \V,  Lafayette,  By  47907 

•*  Mailing  address:  Dr.'R.  D.  Mehta,  Mail  Slop  260-1,  NASA  Ames 
Research  Center,  MolTcIt  Field,  CA  940.?.^ 
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ble  1.  For  a  given 
between  the  two 
y  well. 

cross-wire  probe 
a  fully  automated 
controlled  by  a 
:d  in  T-Z  planes 
probe  oriented  in 
ics  were  averaged 
'  400  samples  per 
nponents  of  mean 
le  Reynolds  stress 


Condition 

t'. 

(m/s) 

dq,, 

(cm) 

(t^l 

Re, 

H 

High-speed  side,  base  case 

15.0 

0.40 

0.05\ 

525 

2.52 

Low-speed  side,  base  case 

9.0 

0.44 

0.06 1\ 

362 

2.24 

High-speed  side,  reversed  case 

l.i.O 

0..^9 

0.0.S4  ' 

t.5.^2 

>  “'9 

Low-speed  side,  reversed  case 

9.0 

0.44 

0.055 

\322 

2.61 

